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Using the Stochastic Adhesion Model (SAM) as a simple toy model for cosmic structure forma¬ 
tion, we study renormalization and the removal of the cutoff dependence from loop integrals in 
perturbative calculations. SAM shares the same symmetry with the full system of continuity-l-Euler 
equations and includes a viscosity term and a stochastic noise term, similar to the effective theories 
recently put forward to model CDM clustering. We show in this context that if the viscosity and 
noise terms are treated as perturbative corrections to the standard eulerian perturbation theory, 
they are necessarily non-local in time. To ensure Galilean Invariance higher order vertices related to 
the viscosity and the noise must then be added and we explicitly show at one-loop that these terms 
act as counter terms for vertex diagrams. The Ward Identities ensure that the non-local-in-time 
theory can be renormalized consistently. Another possibility is to include the viscosity in the lin¬ 
ear propagator, resulting in exponential damping at high wavenumber. The resulting local-in-time 
theory is then renormalizable to one loop, requiring less free parameters for its renormalization. 


I. INTRODUCTION 

Understanding the evolution of density fluctuations under the influence of gravity is a central, but still open issue 
in cosmology. While fluctuations on large scales are small and can be well described using linear perturbation theory, 
on small scales the fluctuations grow large and the Standard Perturbation Theory (SPT) fails, see [T] for a classic 
review. Two possible sources of the failure of SPT have been identified. The first relates to the influence of very 
long-wavelength modes, still within the perturbative regime, on very small scales. After mia first presented their 
resummation scheme, a lot of progress has been made towards a better understanding of how long-wavelength modes 
effect smaller scales [iHin], but also different resummation schemes have been explored mnn]. The effect of long- 
wavelength modes on very short modes is absent for equal time correlators due to Galilean Invariance (GI) [T51116j . 
Nevertheless, it is still possible that long-wavelength modes can strongly affect intermediate modes and one can hope 
that the effect of the former on the latter can be inferred from the effect on short modes by adopting a Galilean 
invariant resummation scheme (see the appendix of [3] for one possible approach). 

The second source of failure, that large density fluctuations on small scales can in principle have a sizable effect 
on the small fluctuations on large scales, must be investigated independently of the importance of long-wavelength 
modes. On small scales not only is the density contrast large but also, and crucially, the single-stream approximation 
fails and hence the fluid approximation fails. Calculating the effect of small on large scales is therefore not possible 
within the framework of SPT or Lagrangian Perturbation Theory (LPT). To address the issue, it has been proposed to 
average over small scales such that one is left with equations for the large scale density contrast 6 and velocity ^ = w 
plus an additional effective stress tensor a which encodes the small scale dynamics [IlHIS]. The Euler equation then 
takes the form 

i9dw-I- (w -I- V$) -I- w • Vw =( 1 ) 

2 L) 1 -b 0 

where D is the linear growth factor, the background cosmology is encoded in 7 = and $ is the rescaled 

gravitational potential, determined by the Poisson equation For a ACDM cosmology we have 7 « 1. 

The effective stress cr is unknown. One possible strategy of treating this term is to measure it directly from N-body 
simulations, in which case cr acts as a source for the perturbations on large scales [191 [20]. Another possible strategy 
is to attempt its parametrization in terms of the velocity and the density contrast as follows 

= j_ + X V X w + ... , (2) 

1 -b 0 D 

where the Vi are effective viscosity parameters and J is a stochastic noise. The ellipses denote terms which are higher 
order in fields or derivatives which are expected to be suppressed compared to the leading terms. This second strategy 
goes under the name of Effective Eield Theory of Large Scale Structure (EFToLSS) and has, in its eulerian formulation, 
drawn a lot of attention in recent years pTM28j . including the possibility to resum long-wavelength modes [5S]. A 
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Lagrangian version has also been formulated in [301131]. The unknown parameters vi are to be fitted to observations 
or simulations. 

The added effective terms not only encode the small scale dynamics, they also ensure that physical quantities, 
for example correlation functions, are independent of the arbitrary smoothing scale. These terms could be treated 
perturbatively, employing a power counting where the linear propagator does not involve any damping due to the 
viscosity terms - see for example (33) for this kind of propagator in the model considered here. The EFToLSS approach 


effectively corresponds to such a power counting applied ot the full set of fluid equations. As we argue below, to ensure 
cutoff independence of correlation functions with such a power counting the r.h.s. of ([^ must be non-local in time. 
Galilean Invariance of this non-local-in-time theory also dictates the inclusion of higher order non-local-in-time terms 
of a certain form on the r.h.s. of ([^. As we will see, these terms ensure cutoff independence and renormalization is 
possible and consistent with GI. In the EFToLSS literature, the physical meaning of this non-locality in time is traced 
to the lack of a fast timescale for short wavelengths. It should be noted however that fitting the Power Spectrum 
obtained by an EFToLSS calculation to a fully non-linear Power Spectrum is possible with a similar accuracy using 
both local and non-local counter terms [32j . 

Another possibility is to include the effective viscosity terms in the linear propagator, see ( |45[ ). This power-counting 
scheme was employed in [33] • The authors of [331 have used a closely related approach by including a local viscosity 
and sound speed in the linear fluid perturbation equations, and find that better agreement with N-body simulations 
can be achieved compared with the usual SPT calculation. We show that in this case the theory is Galillean invariant, 
local in time and one-loop renormalizable. A theory local in time is expected if the effective terms are dominated 
by sufficiently small scales, where, according to the gravitational free-fall time AD ^ the typical time scale 

is much smaller then the time scale on large scales. This is in particular the case if the cutoff of the theory is given 
by the scale where multi-streaming becomes relevant. For example the viscosity is then given by the microscopic 
viscosity plus a contribution from scales beyond the cutoff. 

Let us emphasize that the effective terms aim to encapsulate the influence on large scales of the highly non-linear 
evolution of short wavelength perturbations, accounting for short wavelength deviations from a single stream fluid. 
They are not meant to include the effect of a possible non-trivial background phase space distribution or an initial 
deviation from a single stream fluid, as is the case for free-streaming particles like neutrinos. Nevertheless, the fluid 
description of such free-streaming particles on scales larger then the free-streaming scale feps can also be interpreted 
as an Effective Field Theory with k-ps as the cutoff. Let us briefly sketch this idea. For n > I moments of the velocity 
distribution are of the order iw" ~ and thus can be treated as perturbative corrections to 5 and w. The first 
corrections are then given by the stress tensor, which obeys, neglecting the third moment. 


doO' + + wV • (T 


cr • Vw -f (V’w)’^ • (T = 0. 


(3) 


Splitting the stress into a background and a perturbation u = a\ + 5 (t one can write the stress, for 7 = I, as 


<T=ld.( + rdry(Vw + (Vwn(7) + ..., 


D 


D 


D 


(4) 


' Di 


with Gi and 5(Ji being the background value and the perturbation of the stress at some initial time Di. Plugging this 
into the right hand side of the Euler equation we have 


—TTj—= 


CT, I ^ 1 V.5 ■ 


1(7,; 


D 

-jj] J dr] + W ■ Tiv){r]) + - 


(5) 


To this order the initial stress perturbation plays the role of a stochastic noise and the background stress induces a 
local sound speed and a non-local-in-time viscosity. Higher order contributions can be obtained straightforwardly in 
the double expansion in and the fields S and w. Note that the time dependence of the effective terms is fixed 
and does not coincide with the SPT loop time dependence so the theory cannot be renormalized. Since the velocity 
moments of the background and the initial distribution are in principle known, they can be resumed such that one 
obtains a theory which is non-local and valid at all scales as long as the density contrast is small [33]. We see that 
the effective long wavelength theory for GDM may be thought off as analogous to this approach to free-streaming 
particles. The analogy is imperfect though, given that different time dependences of the effective terms are required 
for GDM. 

The outline of our paper is as follows. In section [ll| we introduce our simplified toy model, the Stochastic Adhesion 
Model. Then in section [lTT| we discuss GI and how it constrains the allowed terms that can be used to parameterize 
the effective stress tensor. In section HVl we discuss renormalization of the local-in-time and non-local-in-time versions 
of the effective theory. We conclude in section [V] 
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II. THE STOCHASTIC ADHESION MODEL 


Instead of discussing the full set of equations consisting of the continuity and the Euler equation, we study the 
technically simpler Stochastic Adhesion Model (SAM) as a toy model. The SAM, as already discussed in [33] . 
can be obtained from the fluid equations by a Zeld’ovich approximation, see also |36H39j for earlier work on the 
Adhesion Model and the Burgers Equation. The Zel’dovich approximation reads w = —and as a result dw/dD = 
9£)W + w • Vw = 0. This approximation decouples the Euler equation from the continuity equation. In principle one 
can obtain 6 from the continuity equation once h is known, for details see [33] . In the following we will not consider the 
continuity equation, since the scope of this paper is to discuss the interplay between renormalization and GI, which can 
be done by considering the Euler equation alone. The SAM is obtained by writing dw/dD = ^^jW+w-Vw = 
and expressing the effective stress tensor as discussed above. It thus approximates deviations of fluid elements from 
their long wavelength Zel’dovich trajectories. One ends up with a time dependent Kadar-Parisi-Zhang (KPZ) equation 
[40] for the velocity potential h, defined by w = — V/i. 

SAM is a stochastic field theory and is most conveniently formulated by defining the MSRJD action^ 


S=l [ dDdD'^P^ 

2 J (27r)3 


. . f 0 {-aD + i'k^)S{D-D') 

(hk, Xk)jj I iM{k,D,D') 




int 


( 6 ) 


1 


IDD' 


d^k 


( ?ik . Xk ) Gq ^ ) + Cint 


where 


and 




A-1 _ f ° [*3^1 

0 “ ^ iM{k) 


^int — 


d^qi d^q2 

(27r)3 (27r)3 


(qi-q2) Xk^q./iq, ^(k + qi+q2) 


(7) 


( 8 ) 


is the interaction vertex. In the second line of Q we condensed the notation for time integrations. The field x is an 
auxiliary field and the term in the action encodes stochasticity. It can be used to encode both stochastic initial 
conditions as well as the continuously acting stochastic part of the tress tensor, modeling the action of small scale 
fluctuations. Accordingly, Af contains the initial Power Spectrum as well as the Power Spectrum A = (JJ) of 
the (gaussian) noise J. For the Fourier transform of Af we assume 


J^{D, D'- k, k') = P$.„(fc)(27r)35^(k + \d)5D{D - D,^)5 d{D' - A„) + A(P, + k'). 


(9) 


with the scale independence of A ensuring that the small scale fluctuations induce the well known peculiar velocity 
Power Spectrum at large scales and, correspondingly the kf^ tail in the density. The term uV^h is the effective 


R{A) 


is used to denote 


viscosity term. G^^^^ is the free Retarded (Advanced) Green function and the notation G, 

the operators appearing in § with Retarded (Advanced) boundary conditions. All correlators of interest can then 
be obtained from the generating functional 


Z = / DhDxe 


tS 


( 10 ) 


with the MSRJD propagator G defined as the functional and matrix inverse of the matrix in the quadratic part of 
the MSRJD action: 


Go = 


(J 


Fo(D,D'-,k) -iG^{D 


-iGf{D,D'-k) 0 


,D'-,k)^ 


( {h^(D)hl(D')) {h^(D)xl{D')) ' 


\{XkiD)hl(D')) {xkiD)xUD')) J ■ 


( 11 ) 


The advanced and retarded Green functions are not independent since Gq {D,D') = Gq*{D',D). Furthermore, 
knowledge of Gq allows the computation of Fq as 


D 


Fo{D,D';k)= J dudv GQ{D,u-,k)Af{u,v,k)GQ{v,D';k). 


( 12 ) 


^ Named after Martin, Siggia, Rose, Jansen and De Dominicis m- 













4 


ki 


G^{D,,D2-,k) = Di 


- D 2 Fo(Di, D2', k) — Di 


Do 


ki ■ ko 


FIG. 1: The Feynman rules for the action § - for more details see [33]. As usual, wave-vector conservation applies and the 
vertex has to be integrated over internal momenta and time. 


As we already discussed, the theory is expected to be non-local in time and this can be implemented by the 
replacement v{D)V^h{D) —>■ dD' v{D,D')V^h{D'). From the action (§ one can read off the Feynman rules 

depicted in figure[^ Note, that the form of the propagators Gq{D, D'\ k) and Fq{D, D'; k) is different if the viscosity 
and the noise are treated pertubatively (a la EFToLSS) or non-perturbatively. Note also that Fq is simply the linear 
Power Spectrum Pl of h. 

The MSRJD action and propagator imply that the ’’self-energy” (the sum of all IPI diagrams^) also has the 
structure 


( 0 

^ i^{D,D'-k) J ■ 


(13) 


If G is the full dressed Green function it satisfies the Schwinger-Dyson equation 

(Go-i-E)oG = i (14) 

where the circle product denotes integrations over time and matrix multiplications. Writing it explicitly we obtain 


{dD + i^k^)G^{D,D')- 

J du E^(G, u)G^(m, D') =5{D-D') , 

(15) 

{-dD + >yk^) G^{D,D')- 

J duY.^{D,u)G^{u,D') = 5 {D - D') , 

(16) 

{do + vk^) F{D, D') - J du E^(n, u)F{u, D') + J 

du {M{D, u) - $(G, u)) G^(u, D')=0. 

(17) 


To close this section, we emphasize that SAM is not intended as a tool for precision cosmology calculations as 
it contains uncontroled approximations. Nontheless, it can reproduce qualitatively the morphological structure of 
the cosmic web obtained from N-Body simulations - see for some early works on the adhesion model. The 

addition of a stochastic component could be used to parametrize short scale, highly non-linear processes and could 
improve the results of those early works by generating more realistic short scale power. Furthermore, it seems that 
an irreducible stochastic component is necessary for describing the effects of short scales and becomes dominant 
over further additions to the stress energy tensor |3S|. With these remarks we postpone a detailed evaluation of 
SAM for future work. What is important for us here is simply that the SAM is invariant under Extended Galilean 
Transformations (GT), like the complete fluid equations, and we are therefore able to discuss the interplay between 
Galilean Invariance (GI), non-locality in time and renormalization with the SAM as a conceptually useful toy model. 


III. GALILEAN INVARIANCE 

Symmetries constrain the allowed terms to be added to the fluid equations for CDM. In the previous section we 
already stressed that the fluid equations as well as the SAM are invariant under GT. The GI of the fluid equations in 
the context of LSS and the corresponding consistency relations were already discussed in Glofthe 

fluid equations is the symmetry of the relativistic equations with a non-trivial Newtonian limit |48) and are therefore 


IPI-diagrams are those diagrams which cannot be cut into two by cutting a single line. 
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related to diffeomorphism invariance of the full relativistic theory. A GT is a time dependent boost with a velocity 
(3{D). The coordinates then transform according to 

D' = D (18) 

X —)• x' = X + f dry /3 (77) = X + T (D ), 

J Di 


and the velocity potential transforms accordingly as 

/i(£),x) ^ h{D\x') = h{D,x + T{D)) -X- j3{D). 
The action ([^ then transform as follows 

S')/*, x] -t S[h, x] + SSj^[h, X]+ J dDd^x x • dof^ + 


(19) 


( 20 ) 


where we used that time derivatives are not Gl, but the convective derivative do + is GI. Compared to the 

action in equation equation ( |20[ ) contains two extra terms. The term SSx'[h,x] will vanish for a Gl stochastic 
noise. The term J dDd^x x (~x • i9d/ 3 + 1/3^) contains two unobservable contributions. The f3^ terms simply adds 
a constant contribution to the velocity potential, while the x • do(3 term is a frame fixing term and ensures that the 
homogeneous mode of the velocity is given by (3 m- 

As discussed above, higher order terms are allowed but GI only allows terms built of second or higher derivatives 
of h. For example we can add a second order viscosity 




( 21 ) 


or new vertices as 


A (V^h)^ + g {d,djh) (d^d^h) . (22) 

Let us now have look at the noise and how it is constrained by GI. The Power Spectrum of a statistically homogeneous 
and isotropic noise J is of the form (J(I3i, Xi) J(Zli, X 2 )) = A(I3i, £> 2 , |xi — X 2 I). The noise term in the action is only 
GI if the noise Power Spectrum is invariant 


A{Di,D 2 , |xi - X 2 I) = A{Di,D 2 , |xi + T{Di) - X 2 - T(D 2 )|), 


(23) 


which is only the case if the noise is temporally white 


A{Di,D 2 , |xi - X 2 I) = A(Di, |xi - X2|)<5 d(7Ai - D 2 ). (24) 

Similar arguments hold for higher order correlators or a multiplicative noise. 

A noise with a finite correlation time is not apparently GI. However, consider the noise term evaluated along the 
path of a fluid element [23] 


J(i3;xfl (A A)), 

where the position of the fluid element can be obtained by solving 


l-D 

Xfl (U, £)') = X + / dry Vh(ry,Xfl(£), ry)). 
JD' 


(25) 


(26) 


Since xg does not change under GTs, neither does the argument of J and its correlators are invariant. It is interesting 
to note that a noise of the form given in equation (251 is nothing but a solution of an equation of the form: 


dJ 

Id 


= doJ + 'w -VJ = ..., 


(27) 


i.e. of a derivative taken along the fluid flow lines. The ellipses denote possible further terms consistent with Galilean 
symmetry. So a non-local noise can be seen as new degree of freedom, governed by (27), which can be added to the 
set of equations. As we will discuss in section]^ the time dependence of the noise is fixed by the time dependence of 
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the loops. This means that it is sufficient to provide initial conditions to specify the noise, suggesting that a non-local 
noise arises from coarse graining over the initial conditions.^ 

If we treat the dependence on the fluid path perturbatively we find for the noise 

J{D-iCfi{D,Dij) = J{D-,x) + [ dr;V/i(7 ;;x)-VJ(£i;x)-f .... (28) 

JDi 

For the corresponding Power Spectrum we find a temporal non-white noise plus corrections in form of a multiplicative 
noise 


A (£>i,T>2; \yi^{Dx,Di) - xfl(i:)2, A)|) ^ 

l>Dx 

+ d ?7 (V/i)(r 7 ;xi) • ViA (£>1,1)2; |xi 

JDi 


A(Di,i:)2;|xi -X 2 I) 


X 2 I) + 



dr] (V/i)(? 7 ; X 2 ) • V 2 A (Di, £> 2 ; |xi 


X2|) + .... 


(29) 


So a non-local noise is allowed as long as it is the first term in a series of terms all with the same coefficient function. 
Evaluating the fields along the fluid path Xfl we can likewise generalize local terms containing h to non-local terms. 
For example a non-local viscosity reads 



dD' v{D, £>')V2/i(£>'; xg (£>', A)) 


oD 



dD' v{D,D')V^h{D';x) 

.D' 

dD' / dr] v{D, A)V2 ((Vh(? 7 ; x)) • Vh(A; x)).... 

•)£>in 


(30) 


Similarly to the noise, a non-local viscosity is allowed as the first term in a series of terms, all with the same coefficient. 

Observe that these extra terms can and must contain the velocity itself, so terms with only one derivative acting 
on h appear in the action. The non-local-in-time terms in (29) and (30), a consequence of GI, lead at one loop to the 
new vertices 


£>1, ki 



O 2 , k 2 


/c^ki • k 2 
2 


{u{D,Di)0{Di 


D2) + u{D,D2)e{D2-Di)) 


(31) 


£>1, ki 

/ 

D, k - i = (k • k2 A(A, A; fc2) + k • kiA(A, A; h)) (0(A - D) + 0(A - D)), 

\ 

\ 

A,k2 (32) 

where as usual “momentum” conservation is implied. In contrast to the usual vertex one has to integrate over all 
three times £>,£>1 and £> 2 . In the EFToLSS approach these terms will only appear at second order, so are relevant 
for the one-loop bispectrum and the two-loop Power Spectrum. 


IV. RENORMALIZATION 

We now discuss the implications of GI for the renormalization of UV-divergences. By a UV-divergence we refer to 
the leading contribution from hard loop momenta, irrespectively of whether the loop integrals are finite or infinite. 
In both cases these are unphysical contributions and must be removed from physical quantities by counter terms 


3 


Our discussion of free-streaming particles in the introduction provides another example where the initial stress tensor has been integrated 
out. 
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corresponding to the effective terms. Since in SPT loop integrals are finite for realistic initial conditions, actual 
divergences are not present at low orders, but will arise at higher orders from loops containing the effective terms. 

When renormalizing the loop integrals one has to pay attention to the fact that they have a non-trivial time 
dependence, so the counter terms must match the time dependence of the UV-divergences, either local or non-local 
in time. Since the time dependence of the effective terms is not constrained by GI or any other symmetry, the 
time dependence can always be chosen as required for the cancellation of UV-divergences, but if the effective terms 
are non-local in time the same divergence appears at higher order again. Since the counter terms appear again at 
higher order, these divergences are automatically renormalized. We will explicitly show that this happens at one loop. 
Afterwards, we will discuss how the Ward-Identities (WI) from GI ensure that this will happen at any loop order. 


A. Non local in time 


The approach taken in the EFToLSS leterature treats the effective viscosity and the noise as perturbative corrections. 
The diagrammatic expansion in the EFToLSS therefore consists of the SPT diagrams plus counter-term diagrams for 
the effective terms. In the framework of SAM this translates into the free propagators reading 


G^{D,D'-k)=e{D- D') 
FoiD,D'-,k)=P^Jk), 

where we used initial conditions in the past at Din 0. At one loop the self-energies are 


(33) 


E^(A,U2;fc) = Ui,k --- 


^iDi,D2;k) = ^k) = D,,k - -- 


D2,-k =-^0{Di-D2)l2 



-- Dn.-k 


= - y V2(A:) + ^ro(fc), 


where we defined 


In = 


d?q 

(27r)3 

dP'q 






q-5 




q -I- — 


(34) 

(35) 

(36) 


As a consequence of causality, inherent in the MSRJD form of the action, E is only non-vanishing for D 2 < Di, 
while $ is always non-vanishing. Since the self-energies are not sharply peaked for any Di and D 2 nor for Di « D 2 , 
they require non-local counter terms. It is immediately clear that the dependence of external wave vector k of the 
leading UV-divergences of E matches the k'^ dependence of the viscosity. The leading UV-divergence in $ comes from 
I 4 , the term with the highest power of q under the integral. Taking the limit where the loop momentum q is much 
larger then the external momentum k, q ^ k, we see that the divergence can be absorbed into the scale-independent 
noise. At higher order it is ensured by the form of the SPT kernels, as dictated by momentum conversation, that the 
/c-dependence of the divergences always matches the one of the viscosity and the noise mm- 
Using the above one-loop self-energies the Schwinger-Dyson equations give 


doG^iD, D') + k^ J du (D, u) + G^ (it, D') = S {D - D') , 

0 

D D 

doFiD, D') Fk"^ J du (^ly {D, u) + F{u, D') + J du u, k) - 4>(k)) G^(u, D') = 0 . 


(37) 

(38) 


0 0 
They can be renormalized by constant contributions to the viscosity and noise 

v{Di,D 2 ) = cl 
A(Gi,U 2 ) = A. 


(39) 











_D,k 


I 





D2 , \^2 


FIG. 2: The 6 one-loop contribution to the bispectrum. The upper three renormalize the vertex. The two diagrams in the 
second line generate a multiplicative noise. The diagram in the last line generates a bispectrum for the noise. 


It has been noticed [261 El] that at lowest order a non-local viscosity can be mimicked by a local one. In principle 
this is still possible at higher order but then new terms are necessary correcting for the error made by using a local 
viscosity instead of a non -loc al one. For the noise a similar procedure is not possible. 

As discussed in section III the same constants and A appear as coefficients of vertices, so the same divergences 
as in the vertices must be present in the vertex corrections, which are given by the sum of the three upper diagrams 
shown in figure 


n(D, Di, D 2 ; ki, k 2 ) = -9(0 - Di)0{Di - D 2 ) 


fc^ki•k 
\2 


" h + 9{D - D^)e{D - D 2 ) I 2 + (Di O £> 2 ). (40) 


Note that contributions to 11 which would renormalize the vertex cancel among the three different diagrams. The 
first term in equation (401 has exactly the from required by GI, the second term can be renormalized by a non-local 
version of the g-term in equation (22) while the non appearance of a term corresponding to A in (22) is a consequence 
of the simple form of the vertex. 

The second class of diagrams we have to consider are those generating a multiplicative noise (the two diagrams in 
the second line of figure]^. Their sum is 

DuD2; ki, k 2 ) = 0(L»i -D)( - - ^^l4(fe) + ^k • (¥ 4 (^: 1 ) + Y4(fe)) • k 


16 


fciki ■ k - (ki ■ k)^ 
32 

^fk-ki-fc?(k-ki)^ 


42(fci) -|- 


16 

fc|k 2 • k - (k 2 • k)^ 


32 


Y2{k2) - -k • {klY2iki) + klY2{k2)) ■ k 


Yoiki) - 


fc|k-k2-A:|(k-k2)^ 


40(^2) ) + {Di o D2), 


128 


128 


(41) 
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where we used that / d^q qq^P,;>.^ (|q ~ 11 ) (I*! + l|) = 0 and defined 


Y„+ 2 (fc) = I 


d^q 

(27r)3 


qqg”P$i, 


Pi. 


q + 


(42) 


The leading divergences are again those of Y 4 and Y 4 , where due to rotational invariance Y must be of the form 


Y„+ 2 (k) = F{k)l + G{k)kk. (43) 

Since in the limit q ^ k, Y„_|_ 2 (k) depends only on the absolute value k, we conclude that the leading divergence 
stems form F{k), which is the same as that of I 4 . So we can simply replace F by —In that limit 4' reads 

4-(P,Pi,P2;ki,k2) ~ e{D, - D) - ^^Y4{k2) + ^ (Y4(fci) + Y4(fc2))) + (A ^ C 2 ). (44) 

We are left with two divergences. First, the one in kF'Y^iki) which can be renormalized by a multiplicative noise of 
the form JV^h and the one in k; • kY 4 {ki) which as required by GI is the same as in the Power Spectrum. 


B. Local in Time 


The power counting of the effective viscosity and noise terms discussed above, corresponding to that employed in 
the EFToLSS literature, implies that loops are renormalizable if the effective terms are non-local in time. We will now 
demonstrate that by including the effective viscosity and the noise in the propagators, loops can also be renormalized 
in a theory local in time. In that case the theory is similar to an ordinary viscous fluid with stochastic noise. 

Let us in the following write the viscosity as !z(_D) = c^v{D), with v{D) = 0(1). We chose the time dependence of 
the noise to be A(I?, D') = /S.v^{D)5d{D — D'). The linear propagator and the linear power spectrum are then given 
by 


0^(0, O'; k) = 


(45) 


Fo(0, O'; k) = {k) + A ^ dV 

Jo 

As already argued in [33] on large scales we recover the usual linear Power Spectrum, while on small scales the Power 
Spectrum is dominated by the noise. 


oo/c7-> c7->/ \ /*min(D,D^) 

Fo(0,0';fc) = e"'=’''= A dr, JJ 

Jo 


_ (/cf dv dri 


2c?,A;2 


^min(D,DM i 

0 dr, 


JS' ^('')6»(0 - O') + {D^D') + 0 ((c^fc)-^) 

2c.,k 


(46) 

(47) 


To obtain the last line we performed a partial integration and neglected terms which are exponentially suppressed. 
Repeated partial integration would allow to calculate the O ((c„fc)“^) terms and higher. Note that on small scales 
the Power Spectrum is exponentially suppressed unless O « O'. ^ 

The leading UV-divergences of the self energies are then 


S(0,0 ';fc)~^ 

a 2 r 

ci>(0,0';fc)~^y 


A f d^q {q‘^ + q-k)ci-k _^^2 


2c2 j (27r)3 g2 

^ d^q [q^-'ir 


v\D')e{D - D')G^{D, O'; g)G«(0, O'; |q + k|) 


.0 ^ 3 , ^k^ 2 . ^,,i>\D')0iD - D')G^{D,D';\q-^\)Gl^{D,D';\q+^\) + (D ^ D'). 

iq+^y{q -2 2 

(48) 


More precisely it is suppressed unless c^/c^ drji>{ri) 1, which translates for monotonic (growing) i^(-D) into {D — D') « ^2 
1 
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The above expressions are non-vanishing only for Di « D 2 and they can therefore be approximated as local in time. 
Integrating over time we obtain to leading order 


E (U, D'; fc) 
a>(D,D';fc) 


- D') 

v^{D2)^5d{D - D') 


d^q -I- q • k)q • k 

(27r)3 g2(g2 + (q + k)^) 

(2^)3 (q+k)2(q_k)2(^2 + ^)' 


(49) 


We immediately see that the time dependence matches the one of the viscosity and the noise.^ The g-integrals are 
actually UV-divergent and must be regularized by an appropriate prescription. From inspecting the integrals we see 
that E has a quadratic and a linear divergence. The quadratic vanishes due to rotational invariance and the linear 
divergence has the correct k dependence to be absorbed into the viscosity. There is also a logarithmic divergence 
stemming form the sub leading terms 0{-^), we dropped this term since the divergence vanishes. The linear divergence 
in $ can be absorbed into the noise. 

Observe that the divergences are the same as in the time independent theory with iz = 1, as are the leading 
divergences at higher orders. There will also be additional divergences with new time dependencies at higher orders 
which will require new counter terms, including terms with additional time derivatives. This is not a problem since 
the theory is non-renormalizable and new counter terms must be added any way. It is interesting to note that if 
instead of a scale independent noise we chose one with A ^ in the UV, the theory becomes renormalizable by 
power counting. 

If we follow the same procedure as for the self-energies for the triangle diagrams we find 


n{D,Di,D2-, ki,k2) ~ -^^d{D - Di)6d{D - D 2 ) I 


-b 


^{D, Di,D2] ki, k2) ~ — 16^5 ^ {D — Di)Sd{D — D 2 ) j 


d^q ki ■ qk 2 ■ (ki -b q)q ■ (k -b q) 

(27r)3 g2(g2 q)^)(g2 -b (k - q)^) 

d^q ki ■ qk2 ■ q(ki-bq)(k2 - q) \ ^ n i 

(2^)3q2(g2 + (k,+q)2)(g2 + (k2_q)2)y 1 2 I 

d^q {kf - 4g^) k • (ki -b 2q) (ki - 2q) • (ki -b 2k - 2q) 
(ki -b 2q)^ (ki - 2q)^ 

^ (50) 


(2^)3 


-b 


-b 


(fc^ -b 4g2)(fc^ -b 4g2 -|- q . (ki — 2q)) 

1 


(2fc2 -b (k - 2q)3 -f k • (ki) - 2q)((fc2 -b 4g2) -|_ q . (kj^ - 2q) 
1 


(fc3 -b 4g3)(2fc2 _|_ (k — 2q)2 -b k • (ki) — 2q) 


+ {Di D 2 ) , 


The individual diagrams contributing to 11 are linearly divergent, but the divergences cancel among the integrals, 
ensuring that the vertex is not renormalized, as required by GI. 4* is by power counting logarithmically divergent, 
but this divergence vanishes due to rotational invariance. 

The third triangle diagram with three external dashed lines is finite and no non-gaussian noise is required to 
obtain finite results. Similarly, higher order correlators are finite at one-loop, so the SAM and most likely also the 
fluid equations are one-loop renormalizable. Of course, the loop integrals still contain unphysical contributions from 
modes beyond a UV-cutoff ki and must in principle be renormalized, but the error made by not renormalizing 
the triangle diagrams is suppressed by ki/A and therefore of the same order as the residual cutoff dependence of 
renormalized self-energies, and hence only subleading. So the minimal set of counter terms required to make the 
local-in-time theory cutoff independent at leading order, and hence the number of free parameters, is smaller than 
in the non-local-in-time version. Note that finite contribution from the g- and A-type vertices or new noise terms as 
discussed in section [ITT] can still be important but a detailed analysis is beyond the scope of our paper and we leave 
it for future work. 


^ Reference m treated the special case of iy{D) = D. 
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C. Ward Identities 


The Ward Identities (WI) encode the statement that the effective action T transforms under infinitesimal GT in 
the same way as the bare action S. This implies a set of relations that counterterms necessarily satisfy. The effective 
action T is related to the generating functional of connected correlation functions and any physical information about 
a system can be obtained from it. We will now briefly discuss how the WI ensure that in a non-local theory the same 
divergences arise in different n-point functions. 

The WI for the fluid equations were already derived in HUES], see also [JS]- The corresponding WI for the SAM 
are 


6T 


ST 


k . T,D) 


+ 5 , 5 (k) ( P{D) ■ 5 k- do^iD) ■ duXuiD) + '-xUD)P\D)'^ 


= 0 . 


(51) 


The 9kX"t6rm relates the mean velocity to the Galilean boost /3 and the y-term is related to an unobservable shift in 
the potential. These terms correspond to the change of the bare action under GT, see equation The 1-particle- 
irreducible vertices (IPI-vertices) defined as 


p(n,m) _ 


jn+mp 


^Xpiim) ■ ■ ■ SXp^MShqi^Xi). ..Shq m ('^m) 


(52) 


/i=X=0 


are given by the sum over all IPI-diagrams with n dashed and m solid lines. Taking n derivatives with respect to y 
and m derivatives with respect to of equation (51), we arrive after a partial integration, for n ^ 1 and m 7 ^ 0, at a 
relation between a vertex with n -I- m legs and one with n -|- m -I- 1 legs of the form 

n m 

Y, V^5D{D - -t- ^ ci,5d{D - 

- J d^kSo (k) 9 i 59 kr (”>™+5 = 0 . ( 53 ) 

Note that = 0 for any m. The WI in this form relate a IPI-vertex to another IPI-vertex with one more leg. 

The extra leg thereby corresponds to a velocity ie a solid line. The delta-function picks out the limit A: —>■ 0. Terms 


with more then one factor of k are therefore unconstrained by the WI, as expected from our discussion in section III 


Since each velocity leg must contain at least one factor of k the WI relate IPI-vertices with velocity legs with one 
factor of fc to a IPI-vertices with one velocity leg less. This is exactly what is required by GI as discussed in section 

imi 

Let us illustrate that this ensures that loops can be renormalized with GI counter terms at any order in perturbation 
theory with the example of 2-vertices and 3-vertices. For n = I and m = I the WI read 


kiG ^{Di,D 2 ;ki) {6{r] - Di) - 0{r] - D 2 )) = lim (Di, 112, r?; ki, p) 

p—^0 

where the full inverse propagator is given by 

i^2; ki) = {Go\Di,D2; fci) - E(i?i, G2, fci)) 


(54) 


(55) 


and similarly the full vertex is composed of the tree level vertex 7 , containing both the usual vertex and a possible 
contribution from the non-local viscosity, and the loop corrections 


r 5 ’ 2 ) {Di,D 2 ,ri;'ki,'k 2 ,p) = -27 (Gi, £>2, ? 7 ; ki, k2, p) - H (Di, £>2,7; ki, k2, p). 


(56) 


Using that in the local theory the tree level propagator is Gg {Di,D 2 ',ki) = dDiSoiDi — £* 2 ) + kliy{Di)SD{Di — D 2 ) 
and the vertex is 7 = 6{r] — Di)S{r] — 02 )^^^ in equation (54) we find that convective derivative is GI. If the viscosity 


is non-local, its contribution does not cancel and the vertex gets an extra contribution, as discussed in section m 
For n = 2 and to = 0 a similar relation holds. At loop level we have 

E (£)i,£> 2 ; fci)ki {0 (7 - £> 1 ) - 0 (rj - D 2 )) = lim SpII (£> 1 , £> 2 , 7; ki, p) 

p—>-0 

$ (£)i,£> 2 ; fci)ki {0 (77 - Di) -0{'n- D 2 )) = lim (£>i, £> 2 , 7 ; ki, p). 

p—>-0 


(57) 
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7/. k rj.k 7/. k 



Di-ki D2.k2 Di-ki " D2.k2 


FIG. 3: The WI in equation (571 at one loop. 


As an example the relations at one loop are depicted in fig. If we now use that in the limit p —>■ 0 the vertices are 
n(ki, p) oc p • kin(fci) and 'I'(ki, p) oc p • ki'I'(fci), we explicitly find the same relation as as we found in section III 


n(Z7i,D2,p;ki,p) = p • kiS(Z7i,D2; A:i) {6{r] - £>i) - 6»(p - £> 2 )) + O (p^) 

T'(£>i,£»2,p;ki,p) = p • ki^{Di,D2;ki) (6»(p - £>i) - 6»(p - D 2 )) + O (p^) . (58) 


These relations guarantee that the same divergence appearing in loop corrections to the power spectrum also appears 
in loop corrections to the bispectrum. In particular since the leading UV-divergences scale as E cx and $ oc 
the WI guarantee that the same divergence also appear in 11 and if the theory is non-local in time, while if the 
divergence is local in time the WI guarantee that there is no vertex correction of this type, unless the self energies 
contain time derivatives. Then as at tree level the WI guarantee that the time derivatives are always part of a 
convective derivative®. 

From the form of the WI in equation (531 it is immediately clear that the same procedure also applies to higher 
order vertices. The general relation between vertices with n -|- m legs and n -|- to -|- i legs can be obtained by applying 
equation (53) recursively. 


D. The Bispectrum in the local and non-local theories 

The scale and time dependence of correlation functions obtained from the local-in-time theory is much more complex 
then that obtained from the non-local-in-time one. But on sufficiently large scales ^ 1, and if the noises are 
neglected, the one-loop Power Spectra of both theories are expected to agree. Nevertheless, the two theories give 
different results for the one-loop bispectrum. In the non-local-in-time power counting, the one-loop bispectrum 
consists of the SPT tree-level and one loop bispectrum plus EFT corrections in form of the new vertex, with the 
coupling g, and the counter terms for E: 


£nLt(£*1, £* 2 , £* 3 ; ki,k2, ks) = £tree(fcl, ^ 2 , fca ) + £l-loop (£^1 , £*2 , D3 ] ki,k2, ^ 3 ) -|- gD ^ 


3 (ki • k 2 ) 


6 


-PL{kl)PL{k2) 


- C„ 


>ki•k 2 


kf 


kn 


PL{k^)PL{k2) -^{Di + SDtDs) + -^{Di + 30^0^) + 


(59) 


permutations. 


® These terms are redundant i.e. for practical calculations they can be simplified using the equation of motion. 
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For the bispectrum in the local theory we find 


fk k ) 

But{Di,D2j Dz]kx,k2,kz) = i?NLT(-Di, Zl2, D^, ki, k2, k^) — gD^ —— PL{ki)PL{k2) (60) 

fc^ki • k2 

+ ^ - PL{ki)PL{k 2 ) + permutations. 

That g appears in i?NLT but not in Buy is an illustration that less free parameters are required for consistency in the 
local theory compared to the non local one. Of course, a g-type vertex can still be added to action but is not needed 
for renormalization. Similiar a A-type vertex can be added to both the local and non-local theory, but is not required 
for renormalization. The second term which is different comes from the extra vertex enforced by GI and by integrals 
in the local theory of the form dr] gf{g) which must be replaced by integrals of the form dg dg' f{g') in the 

non-local theory. These integrals are the same only if / = 1, while for a power law / oc iA" both integrals become 
oc but with a different prefactor. Since for the full system the time dependence of the propagator is more 

complicated one should expect that the difference between Buy and AAnlTi the bispectra in the local and non-local 
theories, becomes larger within the full theory. 


V. CONCLUSION 

In this paper we discussed renormalization of the stochastic adhesion model (SAM). This is a toy model of structure 
formation, based on a simple parameterization of deviations from the Zel’dovich trajectories, which shares the same 
symmetries as the full set of Euler and continuity equations along with extra effective “viscosity” terms and a stochastic 
noise term. Because of these features we expect that the general conclusions arrived at here should also be applicable 
to the complete theory. 

Treating the viscosity terms as counter terms, one is led through a one-loop calculation to a theory necessarily 
non-local in time. To ensure Galilean Invariance, non-local counter-terms must be evaluated along the fluid element 
path, introducing vertices at higher orders which are related to the lowest order counter terms by having the same 
coefficients. These terms are individually not GI, so their dependence on the wave-vector differs from that of local 
counter terms. As we have shown, these terms are required to ensure cutoff independence. This implies that the 
EFToLSS type approach with non-local-in-time counter terms is consistent with respect to renormalization, while the 
EFToLSS with local counter terms is not. 

Alternatively, if the effective viscosity is included in the computation of the linear propagator [3311311 > leading to 
exponential damping on short scales, one obtains a consistent theory, renormalizable to one loop. For a numerically 
small viscosity one naturally recovers the known results with perturbative local counter terms for the one-loop Power 
Spectrum. In contrast to the EFToLSS type of approach though, no new counter terms are required to renormalize 
the one-loop bispectrum in this case. 

We briefly discussed the possibility to decide whether a theory local or non-local in time would be appropriate as 
an effective description of GDM by considering the bispectrum. As we pointed out, the precise size of corrections 
stemming from virialized scales are crucial in that respect. The local-in-time approach requires a minimal set of 
counter terms for renormalization and as a consequence the one-loop bispectrum is fully predicted by the one-loop 
Power Spectrum, although extra parameters are allowed. The non-local-in-time approach requires one more parameter 
for the renormalized bispectrum. We expect that similar conclusions would hold in the full system of continuity -I- 
Euler equations. Renormalized bispectra have only been computed in the EFToLSS approach which uses the non- 
local-in-time power counting. It would therefore be interesting to contrast the predictions of the local and non-local 
in time approaches in the full theory against the bispectrum from N-Body simulations. 
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